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ABSTRACT 


We present the results of a search for astrophysical sources of brief transient 
neutrino emission using IceCube and DeepCore data acquired between May 15th 
2012 and April 30th 2013. While the search methods employed in this analysis 
are similar to those used in previous IceCube point source searches, the data 
set being examined consists of a sample of predominantly sub-TeV muon neu¬ 
trinos from the Northern Sky (-5° < 6 < 90°) obtained through a novel event 
selection method. This search represents a hrst attempt by IceCube to identify 
astrophysical neutrino sources in this relatively unexplored energy range. The 
reconstructed direction and time of arrival of neutrino events is used to search for 
any signihcant self-correlation in the dataset. The data revealed no signihcant 
source of transient neutrino emission. This result has been used to construct 
limits at timescales ranging from roughly 1 s to 10 days for generic soft-spectra 
transients. We also present limits on a specihc model of neutrino emission from 
soft jets in core-collapse supernovae. 

Subject headings: neutrino astronomy, neutrinos, GRB, supernova, astroparticle 
physics 



1. Introduction 


The nascent field of high-energy neutrino astronomy opens the possibility of answering 
several open questions in astrophysics due in large part to the neutrino’s ability to escape 
the densest regions of astrophysical environments. Specifically, the detection of transient 
astrophysical neutrino sources will help shed light on the acceleration mechanisms at work 
in some of the most energetic phenomena in the Universe such as gamma-ray bursts, 
supernovae, and active galactic nuclei. Previous attempts to detect such sources with the 


IceCube Neutrino Observatory (Achterberg et ah 2006) are most sensitive to neutrino 


fluxes above 1 TeV with poor sensitivity below 100 GeV. Searches for astrophysical 


sources at lower energies (1-100 GeV) have been performed by Super-Kamiokande (Thrane 


et ah 2009), however the detector’s 50 kton instrumented volume limits its sensitivity 


to astrophysical neutrino fluxes. A newly developed 30-300 GeV muon neutrino sample 


collected by IceGube and its low energy extension DeepGore (Abbasi et ah 2012b) enhances 


IceGube’s sensitivity in this under-explored energy range. In this paper we will present the 
results of a search for transient neutrino emission in this GeV-scale neutrino sample. 

The detection of astrophysical neutrino sources is a primary design goal of the IceGube 


Neutrino Observatory (Achterberg et ah 2006). Located at the geographic South Pole, 


IceGube utilizes the clear Antarctic glacial ice ice cap as a detection medium for the 
Gherenkov light produced by secondary products of neutrino interactions. The detector 
consists of 5,160 Digital Optical Modules (DOMs) distributed among 86 cables or “strings” 
to form a 1 km^ instrumented volume. These DOMs house photomultiplier tubes (PMTs), 
to detect Gherenkov photons, as well as digitizing electronics for initial processing of the 


PMT data (Abbasi et ah 2009). A centrally located region of denser instrumentation 


featuring DOMs with more sensitive PMTs comprises the DeepGore sub-array. This 
extension to the IceGube array enhances the detector’s response to lower energy neutrino 
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events. 

Typical searches for astrophysical sources with IceCube make use of a sample primarily 
comprised of an irreducible background of high-energy atmospheric muon neutrinos {E^ > 1 


TeV) to look for both steady (Aartsen et ah 2014b) and transient sources (Aartsen et ah 


2015). As of yet, these searches have not found any signihcant self-correlations within the 


data sample nor correlations between the neutrino data and known astrophysical objects of 
interest. So far, these analyses have largely eschewed low energy neutrino events collected 
by DeepCore for two reasons. First, the poorer angular resolution of these events renders 
them less suitable for pointing analyses. Second, the soft spectrum of the atmospheric 
neutrino flux results in higher rate of background neutrino events. However, the increased 
background can be somewhat mitigated by searching solely for transient sources. Therefore, 


applying previously developed search techniques (Braun et ah 2010) to a sample of low 


energy (30 GeV < Ey < 300 GeV) muon neutrino events from DeepGore can enhance 
IceGube’s sensitivity to short transient neutrino sources with softer spectra. 

Due to the large atmospheric neutrino background in this energy range, searches using 
a data set composed of these low energy events will only be sensitive to emission timescales 
on the order of one day or shorter. Active galactic nuclei (AGN) undergoing flaring events 
are one potential source for emission on this timescale. Protons may be accelerated in 
relativistic jets, powered by accretion onto the AGN, resulting in the production of pions 
(and subsequently neutrinos) in shocks due to proton-photon interactions and proton 


self-collisions (Becker & Biermann 2009). For some of the timescales under consideration in 


this search, AGN-powered hadron acceleration must occur over a compact region and will 


require very large acceleration gradients (Klein et ah 2013). The presence of these large 


gradients will result in signihcant acceleration of muons prior to decay, leading to spectral 
hardening of the neutrino hux. Thus, if neutrino emission is occurring over short timescales. 
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it will feature enhanced visibility at higher energies. 

Sub-photospheric neutrino emission from gamma-ray bursts (GRBs) represents another 
possible source for this search. A model for photospheric gamma-ray emission in GRBs 


by Murase et al. (2013) suggests that a substantial flux of 100 GeV-scale neutrinos may 


be produced during the initial stages of relativistic outflow in the GRB. Decoupling of 
protons and neutrons during the initial formation of the relativistic jet causes hadronuclear 
collisions resulting in the production of pions and the production of neutrinos via pion 
decay. The predicted energy for the neutrinos produced in these sub-photospheric collisions 
is on the order of 100 GeV, and therefore this GRB neutrino flux may only be visible to 
IceGube searches with the inclusion of sub-TeV neutrino events. 

Perhaps the most promising potential source for this study is a special class of 


core-collapse supernova referred to as choked GRBs (Meszaros & Waxman 2001). The 


standard GRB model assumes that relativistic jets are generated during the accretion of 


material onto the compact object formed during core-collapse (Rees & Meszaros 1992) 


Fermi-acceleration of charged particles occurs within the internal shocks of these jets leading 
to gamma ray emission once the jets breach the surrounding stellar envelope. There is an 
observed correlation between long duration GRBs and core-collapse supernovae (GG SNe) 


((Woosley & Bloom 2006), (Modjaz 2011)). While the observed fraction of SNe resulting in 


the occurrence of a GRB is quite low, it may be that a larger fraction of core-collapse SNe 
still manage to produce mildly relativistic jets. Due to insufficient energy, these jets fail to 
break through the stellar envelope and any gamma ray emission is effectively ‘choked’ off. 
If protons are accelerated in these jets, then neutrino production will occur in the shocks of 
the jet irrespective of whether or not the jet successfully escapes. A model of this neutrino 


emission proposed by Razzaque et al. (2004) and extended upon by Ando & Beacom 


(2005), hereafter referred to as the RMW/AB model, suggests that these neutrinos may be 
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detectable by IceCube-DeepCore for nearby supernovae (Taboada 2010). Previous IceCube 
analyses have investigated the RMW/AB emission model with respect to a specihc source 


’Abbasi et al. 

2011 

) and as part of the optical followup program ( 

Abbasi et al. 

2012a), 


however the presented search marks the hrst use of low-energy muon neutrino events in 
constraining this model. 


We present the results of a search for transient neutrino emission with a set of 
low-energy neutrino event data collected from May 15th, 2012 to April 30th, 2013. The 
data selection methods used to acquire this unique event sample will be detailed in Sec. 2. 
Analysis methods and search techniques are discussed in Sec. 3. Finally, the results of the 
search are given in Sec. 4 in addition to how these results may be interpreted within the 
context of generic neutrino flares as well as choked GRBs under the RMW/AB model. 


2. Event Selection 

The data acquisition process begins with the fulhllment of one of three trigger conditions 
that prompt the readout of the detector data. Each of these triggers requires some number 
of DOMs to exhibit hard local coincidence (HLC) within a dehned time window. To satisfy 
the HLC condition, two or more neighboring (or next-to-nearest-neighboring) DOMs on the 
same string must register photon hits within a ±1 /is window. The trigger for the lowest 
energy events (often referred to as simple majority trigger 3 or SMT3) requires three HLC 
DOM hits within a time window of 2.5 /is among the DeepCore string DOMs (or in DOMs 
on IceCube strings neighboring DeepCore). The two other triggers that serve as input for 
this event selection operate over the entire detector array with one requiring eight HLC 
DOM hits in a 5 /rs window (SMT8) and the other requiring four HLC DOM hits within a 
cylinder of height of 75m and a radius of 175m in a 1 /is window (Cylinder Trigger). 
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Events satisfying these trigger conditions are then passed to the DeepCore data hlter 


(Abbasi et ah 2012b). This filter reduces the number of cosmic ray muons by using the 


outer regions of the detector as an active veto to tag down-going events originating outside 
the detector. Specifically, the filter examines timing and position information of DOM hits 
inside the DeepCore fiducial volume to identify a center of gravity (CoG) or vertex. For 
each DOM hit in the veto region, the speed of a hypothetical particle connecting that veto 
region hit to the CoG inside the hducial volume is calculated. Veto regions hits whose 
speed lies within a range consistent with that of the speed of light are causally related and 
are therefore likely the product of background cosmic ray muons. Events having more than 
one correlated veto region hit are removed by the filter. 

During the observation period of this search, the DeepCore filter consisted of two 
separate branches characterized by differing dehnitions of fiducial and veto volumes as 



makes use of some isolated DOM hit information instead of only using HLC hits. Events 
satisfying the SMT3 trigger feed the standard DeepCore filter branch whose fiducial and 


veto region definitions are roughly equivalent to those described in Abbasi et ah (2012b). 


The SMT8 and Cylinder Trigger events, in addition to SMT3 events that fail the standard 
hlter branch, feed into the other branch of the hlter which makes use of a more relaxed 
veto region, consisting of two instead of three layers of IceCube strings, providing a larger 
detection volume. The output of both branches of this hlter are used in this search with 
the standard three-layer veto focusing on low-energy events and the two-layer veto branch 
retaining higher energy events. These branches are referred to as the low-energy stream 
(LES) and high-energy stream (HES) and have an exclusive event rate of 17.3 Hz and 23.3 
Hz, respectively. 
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2.1. Veto Cuts and Event Reconstruction 


Events belonging to both the LES and HES are snbjected to several cnts that make 
use of veto region hit information, event topology, and event reconstructions to reduce 
the volume of cosmic ray background events as well as eliminate events that are the 
result of PMT dark noise-induced triggering. The first of these cuts requires at least 
two DeepCore DOM hits within a 250 ns window to remove SMT3 events that are the 
result of spurious hits. An algorithm designed to search for track-like events is then 
used to eliminate noise-induced events that show little evidence of correlation in DOM 
hits. Additionally, events are required to have at least 10 hit DOMs, to allow for a 
well-constrained reconstruction. The DeepCore hlter algorithm is also reapplied several 
times using looser DOM hit cleaning settings to allow more isolated DOM hits in the veto 
region to contribute to the vertex calculations. Finally, the number of DOM hits that occur 
prior to the first hit inside the DeepCore detection volume is used as a cut parameter to 
eliminate potential cosmic ray muon events missed by the filter. 


The initial event reconstruction uses a simple linear £t (Aartsen et ah 2014a) to 
determine the first-guess direction of a muon track that describes the observed DOM 
hit pattern. This linear reconstruction is then used as a seed for a likelihood-based 


reconstruction (Ahrens et ah (2004)) which uses a single-photoelectron (SPE) hypothesis 
to describe the probability of DOMs receiving light from the track at a given time due to 
scattering in the ice. Six iterations of the SPE likelihood reconstruction are performed to 
obtain a best-fit track for the event. Any event with a reconstructed direction, from either 
the linear or SPE likelihood fit, more than 5° above the horizon is removed from the sample. 
We also require that the angular separation between these two reconstructions is less than 
30° for events in the HES sample. 


Spurious DOM hits that occur in the central detector prior to the arrival of cosmic 
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ray muons allow many background events to elude detection through the standard veto 
technique. To isolate these events, a separate SPE likelihood reconstruction is performed 
without using any information from the first two DOM hits in the event. Just as before, 
events with a reconstructed direction more than 5° above the horizon are removed. Events 
in the LES portion of the sample are disproportionately affected by noise hits due to both 
the lower light yield of these events as well as the increased noise rate of the higher quantum 
efficiency DeepCore DOMs. An additional SPE likelihood reconstruction is performed for 
LES events that attempts to mitigate the noise contribution to the likelihood by requiring 
isolated DOM hits to be more strongly correlated to hits satisfying the HLC condition. 
Once again, if the best-fit direction from this additional reconstruction on LES events is 5° 
above the horizon, the event is removed. 


A hnal event reconstruction uses the previously mentioned six iteration SPE likelihood 
fit as its seed. This reconstruction differs from the seed in two important ways. First, it 
uses a multi-photoelectron likelihood (MPE) instead of the simpler SPE algorithm used 


previously (see Ahrens et ah (2004)). Second, a parameterization of Monte Carlo simulation 
of photon transport is used in place of an analytic approximation to model the timing 


distribution for the arrival of Cherenkov photons to the DOM PMTs (Whitehorn et ah 


2013). This reconstruction is identical to that used in a multi-year point source search 


with IceCube (Aartsen et ah 2014b) and the results of this ht are used for the hnal data 
analysis. In order to estimate the angular uncertainty of the reconstruction, the likelihood 
space about the reconstructed direction is ht with a paraboloid via the method described 
in 


Neunhoher (2006). The angular uncertainty derived from the paraboloid method serves 


as event quality parameter, and only events having an estimated angular error ai less than 
45° are kept. 
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2.2. Boosted Decision Tree 


After the application of the described veto and reconstruction cuts, the ability to 
separate the muon background from potential neutrino signal events via simple cuts is 


drastically reduced. We therefore use a boosted decision tree or BDT (Hastie et ah 2001) 


in order to isolate a final sample with acceptable neutrino purity, i.e. < 10% of events are 
the result of background cosmic ray muons. The use of a BDT allows for the classification 
of events by examining several event parameters, and it is a technique that has proven 


useful in previous IceCube analyses (Aartsen et al. 2013). The decision tree is formed 


through an iterative process in which a series of cuts on event parameters is chosen to 
maximize separation between signal and background training samples. Any events that are 
misidentified by the decision tree are then re-weighted or ‘boosted’ to increase the likelihood 
of correct classification by the next iteration of the decision tree. 

At this level of event selection, the large majority of experimental data still consists of 
background cosmic ray muons allowing the actual data to serve as a background training 


sample for the BDT. Simulated muon neutrino events, generated by GENIE (Andreopoulos 


et al. 2010) and Neutrino Generator (a modihed version of ANIS (Gazizov & Kowalski 


2005) specific to IceGube), are used for signal training of the BDT. Neutrino signal events 
belonging to the LES or HES branches exhibit significant differences in the distribution of 
the input BDT parameters, described below, necessitating the construction of two separate 
BDTs. 


The event parameters used for the LES tree include the location of the reconstructed 
event vertex, the number of ‘direct’ DOM hits (featuring a photon travel time residual 
between -25 and 150 ns with respect to the reconstructed muon track), the reduced 
log-likelihood of the MPE reconstruction, the average distance between DOM hits and the 
reconstructed track weighted by DOM PMT charge, and the highest clustering of veto 




















16 


region PMT charge (found by brute force reconstruction methods). The HES BDT makes 
use of the direct hits parameter described above, the reduced log-likelihood of the MPE 
reconstruction, the average charge-weighted DOM distance to track, and the best fit track 
length using information from direct DOM hits. A simulated signal neutrino event sample 
weighted to a E“^'® (LES) or E“^ (HES) spectrum is used for signal training. 




Fig. 1.— Final event rate distributions for the number of DOMs registering hits during the 
event (left) and the cosine of the reconstructed event zenith in detector coordinates (right). 
The solid line describes the final level dataset while the black points represent the sum of 
the various simulated species. Bins featuring large error are the result of atmospheric muon 


events, generated by CORSIKA (Heck et al. 1998), which suffer from limited statistics at 
the final level. 


Events are then input to the trained BDT, and a cut on the event BDT score is 
imposed to yield a data sample featuring a neutrino purity of approximately 90%. This 
final event sample consists of 22,040 events over a livetime of ~330 days, corresponding to a 
data rate of about 0.77 mHz. As Figure [^indicates, the final sample is mostly composed of 
atmospheric neutrinos with an estimated cosmic ray muon contamination of approximately 
0.07 mHz. There is a disagreement between simulation predictions and experimental data in 
the rate of events featuring a low number of DOM hits. The source of this discrepancy is not 
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fully understood, however, the distributions of other event parameters, e.g. reconstructed 
zenith shown in Figure are well-described by the atmospheric neutrino simulation. Given 
the agreement between simulation and data for event parameters relevant to the analysis 
method, we contend that the simulation of signal events for the purpose of calculating 
the sensitivity of the search is accurate. Additionally, we do not rely on simulation for 
modeling of analysis background, and we instead use the experimental data itself to directly 
determine the background characteristics. 

The neutrino effective area for this event selection is shown in Figure]^ While standard 
IceCube analyses clearly have superior sensitivity at higher energies, this event selection 
shows increased acceptance for events below about 100 GeV in neutrino energy. Figure 
also shows the angular resolution for events at the analysis level as a function of energy. 
Lower neutrino energies result in muon tracks that are both shorter and dimmer, leading to 
difficulty in resolving the direction of the neutrino primary. The kinematic angle between 
the neutrino primary and muon secondary also contributes to the angular error. The 
median kinematic muon-neutrino angle after event selection ranges from ~3° at 50 GeV to 
~1° at 300 GeV. As Figure shows, the efficacy of the reconstruction method used in this 
analysis begins to deteriorate rapidly below 30 GeV due to too few DOM hits. Although 
the pointing ability of these low-energy neutrino events is limited, they are still able to 
contribute to the search through temporal correlation with other events in the sample. 


3. Analysis Method 

The search methods employed in the analysis of this data are nearly identical to those 


used in previous time-dependent IceGube analyses (see Braun et ah (2008) and Aartsen 


et ah (2015)). The arrival times and directions of events within the dataset are input to a 


likelihood function which is then used to perform a likelihood ratio test to compare a signal 
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Fig. 2.— (left) The muon neutrino effective area as a function of neutrino energy for the 


presented search. The effective areas for both the 4 year IceCube point source search (Aartsen 


et ah 

2014b 

) and the 4 year ANTARES point source search ( 

Adrian-Martinez et ah 

2012) 


are plotted as well for comparison, (right) Muon neutrino angular resolution as a function 
of energy after event selection. 


plus background hypothesis for the data to the background only hypothesis. 

Construction of this likelihood function begins with the assignment of individual 
event probabilities that reflect the likelihood of seeing an event i with arrival time ti, 
reconstructed direction Xj, and angular uncertainty Uj given a hypothetical source located 
at Xs with strength Us having a Gaussian time prohle with mean time to and width a^- 

'Pi(xi, U, (Ti|xs, Us, to, = —Si + ( 1 - — ) Si (1) 

^tot V ^tot / 

The Si and Bi terms listed in Eq. are the signal and background probability density 
functions (p.d.f.) respectively. The p.d.f.s used in this search differ slightly from those in 
previously reported searches in that they use no reconstructed energy information. The 
signal p.d.f. is given by 


Si (|Xj X^ I, f j, to, (J.UJ, (Jj) Si (|Xj X^ 15 *^*) ' 'Bi (ti, to, , 


( 2 ) 
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where 


and 






K, 


dvr sinh n 


exp (k cos |xj — Xs|) 


■y/^i 


exp 


TTCTqi 


{tj - tpf 

2al 


( 3 ) 


( 4 ) 


The spatial component of the signal p.d.f., S',, is the Kent-Fisher distribution (Kent 1982), 
and it represents a slight deviation in the signal p.d.f. definition with respect to previous 
searches (see Aartsen et ah ( 2014b[ )). This function is analogous to a 2-dimensional 
Gaussian distribution, but it is normalized to the 2-sphere rather than an infinite plane. 
The concentration parameter k is determined by the event angular uncertainty and is 
defined as k = The temporal component of the signal p.d.f., Tj, is simply a Gaussian 


with mean emission time of to and a width of cr^. 


The background p.d.f., i3j, is derived from the hnal level data set which is dominated 
by background. It has the following form 

i3i(Xi,ti) = (5) 

where T is the total livetime of the search, PshgOed^i) is a p.d.f. describing the event 
declination distribution, and PshgAzio^i) is a p.d.f. describing the event distribution in 
detector azimuth. These p.d.f.s are generated directly from data, without reference to 
background simulations. 

The likelihood function itself is simply the product sum of all individual event 
probabilities: 

T (x^, TLg , tg, Pi (\) diPoy ^wi (6) 

The ratio between the likelihood function values under the background only hypothesis 
{ris = 0) and the signal plus background hypothesis is maximized by varying the source 
parameters n^, a^, and to- The test statistic A is then defined as the maximum value of the 
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likelihood ratio: 


A = —2 log 


■\/^aw C,{ns = 0 ) 


(7) 


T 72^, to, 

with Cijig = 0) corresponding to the likelihood of the nnh hypothesis and C{'Xs,ns,io,d'w) 
the likelihood of the signal plus background hypothesis with the best-ht values of the source 
parameters. Because this is a search for sources of hnite duration over a limited timescale, 
the number of potential short duration flares within the data set exceeds that of flares of 
longer duration, leading to an effective trials factor. This results in a bias towards flares of 
shorter duration. We counteract this effect by introducing a marginalization term T/a/^o\u 
in the test statistic formulation which serves to penalize flares of shorter duration. This 
term also ensures that the test statistic will asymptotically follow a distribution with 
degrees of freedom corresponding to the number of htted parameters for data consisting 
solely of background events. More details about this term and its justihcation can be found 
in 


Braun et al. (2010) 


The behavior of the test statistic enables the maximized value A to be used to 
estimate the pre-trials p-value of the best-ht hare through the invocation of Wilks’s theorem 


(Wilks 1938). Because this search attempts to maximize the signal hypothesis over the 


whole Northern sky many times, the actual signihcance of a given hare must be adjusted to 
account for the ehective number of trials accrued during the sky scan. We use the procedure 


detailed in Aartsen et al. (2015) that involves scrambling the event arrival times in the hnal 


dataset, which also serves to scramble the event right ascension. The search is performed 
on the randomized background data set and the p-value of the most signihcant hare in 
the search is recorded. Many iterations are performed to build a distribution of p-values 
which can then be compared to the p-value of the result from the real data. The fraction of 
background trials that result in a p-value of equal or greater signihcance than the observed 
p-value dictates the probability that the observed result is simply the consequence of a 
random background huctuation. This probability is referred to as the post-trials p-value and 
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it represents the true significance of the search result with proper trials factor correction. 

In order to preserve generality, the presented search makes no use of information 
outside of the data set to designate source regions or time periods of interest. Instead, 
each point in the sky over a declination band ranging from -5° to 90° is examined. This 
is accomplished by discretizing the sky into separate bins and letting the location of these 
bins serve as a grid over which to test a hypothetical flaring source. As this is an unbinned 
likelihood analysis, the data itself is not binned and events may contribute to the likelihood 
at any location being tested. Maximization of the likelihood is then performed to obtain a 
test statistic A for each grid point. The first iteration of this scan uses a relatively coarse 
2° by 2° binning. After this first scan, a followup scan with finer 0.5° by 0.5° binning 
is performed over the coarse bins featuring a pre-trials p-value more significant than a 
predefined threshold (—logio(p-value) > 1-75). The result is a map of pre-trials p-values 
which shows the estimated significance of the best-fit flare hypothesis at each grid point in 
the scan. The best-fit flare from the point featuring the most significant maximized test 
statistic after both scans is returned as the hottest spot in the search. 


4. Results and Interpretations 

Applying the described analysis method to the unscrambled dataset yields the skymap 
of the pre-trials p-values shown in Figure]^ The most significant flare is located at (RA, 
Dec.) = (268.75°, 54.25°) with a signal strength of 13.53 signal events and a width 
of 5.89 days with the peak occurring on MJD 56107 (2012 June 29). The pre-trials p-value 
for this flare is estimated at 6.68x10“®. The test statistic for this flare is compared to a 
background test statistic distribution constructed from 2 x 10^ scrambled trials in Figure]^ 
The background distribution gives a post-trials probability of seeing such a flare in a data 
set consisting of background only is 56.7%, indicating that this flare is entirely consistent 
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with the background hypothesis of the data. In light of this null result, we can set an 

+ 90 ° 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 


-logic P 

Fig. 3.— Sky map of pre-trials p-values for best £t flares per bin. The black circle identifies 
the location of the most significant flare found at RA = 268.75° and Declination = 54.25°. 
The nominal resolution of the map is 2°, however regions with more significant p-values 
receive a finer 0.5° resolution. 

upper limit on the time-integrated neutrino flux of any possible unobserved neutrino flare 
that may have occurred during the search period. 


4.1. Generic Source Limit 

Due to the focus on low-energy events in this search, we choose to examine the 
limit with respect to a soft-spectrum E~^ generic flaring neutrino source with a Gaussian 
emission prohle. An upper limit is established through signal injections at a specified 
location through the following process. First, we select the p-value of the most significant 
flare found in the data to serve as a threshold for signal injection trials. Signal events are 
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Background Trials Test Statistic Distribution 



logio(A) 


Fig. 4.— Distribution of maximized test statistic A for 2 x 10^ searches performed on 
randomized datasets. The dashed line indicates the value of A for the most signihcant flare 
found in the data. 

then injected with some Poisson mean value that is increased until the recovered p-values 
from the injections exceed the threshold p-value 90% of the time. This Poisson mean 
number of signal events is then taken as an event upper limit for the analysis method. 

The upper limit on a generic flaring source for several emission timescales and choices of 
declination is plotted in Figure The number of events required rises at longer timescales 
as the rate of accidental background correlations becomes non-negligible. The limit in terms 
of time-integrated flux (GeV“^- cm“^) is also plotted. This limit is obtained by folding the 
source spectrum with the effective area of the event selection and normalizing the flux so 
that the number of events produced in the detector corresponds to the calculated Poisson 
mean event upper limit. 
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Time-Integrated Flux Upper Limit E ^ 




LogiQ{(7^) (Days) 


Fig. 5.— (left) Upper limit (90% C.L.) for a generic transient sonrce as a fnnction of 
flare width cx^ averaged over the declination range of the analysis. The limit is given in 
mean number of events (left axis) as well as in time-integrated flux at a reference energy of 
100 GeV (right axis), (right) Upper limit (90% C.L.) for a generic transient source as a 
function of flare width for different values of source declination. 

4.2. Choked GRB Limits 

This null result can also be used to construct limits on specihc neutrino emission 
models such as the RMW/AB model for choked GRB emission mentioned previously. 
Unlike the hard spectra sources (e.g, E“^) that are the typical target in IceCube searches, 
the neutrino flux for choked GRBs is predicted to be much softer. The spectral shape can 
be modeled via a doubly broken power law with spectral breaks occurring as hadronic 
{E^(i)) and radiative {Ej^( 2 )) cooling mechanisms become efficient (see Eq. |^. Using the 
canonical RMW/AB model parameters, the break energies for pions (kaons) occur at 30 
GeV (200 GeV) and 100 GeV (20 TeV). Therefore the neutrino spectrum is predicted to be 
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very soft at > 1 TeV energies. 


~dE 
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m 


27rDnn{E'p„^^^/ 


( 8 ) 

(9) 


The fluence E^ at the Earth is given by Eq. [^and depends npon the pion (kaon) mnltiplicity 
< n >, the nentrino production branching ratio for pions (kaons) the minimum and 

maximum proton energies the kinetic energy of the jet Ej, the bulk Lorentz 

factor Tf), and lastly the distance to the source D. Equations and reveal that the 
normalization of the neutrino flux at the Earth is highly dependent on the kinetic energy of 
the jet Ej and the bulk lorentz factor Tb. These two parameters also determine the shape 
of the spectrum as the hadronic (Ej,(i) oc Ej^Tl) and radiative {E^( 2 ) oc Tfe) break energies 
depend upon these jet properties as well. We therefore choose to examine the predicted 
neutrino fluence in Ej-Tb phase space. 


To determine which values of these parameters produce a fluence detectable through 
our search method, an event upper limit is first determined via the injection of signal events 
following a spectrum set by the value of Ej and T b (the same process used to calculate event 
sensitivity for the generic E~^ scenario). The emission timescale for these injections is given 
by a Gaussian with a width of 100 s. The calculated event upper limit is then combined 
with the effective area of the event selection to determine the neutrino fluence necessary for 
detection. For a given choice of Ej and Tfe this sets a limit on the distance at which the 
source would still be visible to the search, and we define this distance D.uis as the visibility 
distance. 


When combined with the area of sky examined by the search this visibility distance, 
in turn, dehnes a parameter dependent volume Va (= over which the search 
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method monitors. This monitored volume corresponds to the region in which a choked GRB 
event should be visible to the presented search method with 90% confidence (assuming jet 
alignment). If the observation period of the search is considered, this monitored volume can 
be converted into a limit on the volumetric rate of choked GRB events as a function of Ej 
and Tfo. This, however, requires the assumption that the jets of any choked GRB event in 
this volume are aligned with the Earth. To obtain a limit more representative of the actual 
distribution of choked GRB orientations, one can include a geometrical correction factor 
that takes into account the opening angle of the jets, dj, which is often approximated as 


6j ~ (Mizuta & loka 2013). Because the physics that determines this opening angle is 
not entirely known, we choose not to include any correction for jet opening angle. The rate 
limit is then given by 


V ' 


( 10 ) 


where r is the livetime of the search, Va is the monitored volume previously defined, and 
[/.L.(0j/u) is the null observation upper limit on the number of choked GRBs that occurred 
in our monitored volume with background expectation of /r. 


We define this background expectation /x as the expected number of ‘false positive’ 
flares that occur due to coincident background events during a given search. To calculate 
the value of ju, we first perform many iterations of the analysis with signal events injected 
at a specific declination where is the calculated event sensitivity at that declination. 
The test statistic for these injection trials form a distribution from which we can take the 
median value. A”)®'’*. Once A”)®'’* has been determined, the analysis is run again scanning 
over the same declination band using a time-scrambled dataset with no injections. Several 
iterations of this procedure builds a background test statistic distribution. The number of 
entries in the background distribution whose test statistic value exceeds the threshold A”)®'’* 
are then recorded. This number is then divided by the number of background-only analysis 
iterations performed to yield an expected false positive rate per search. This procedure 
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revealed the false positive rate to be very small (< 10 ^). We therefore take /r ~ 0 leading 
to a Neyman upper limit of 2.3 from the null observation. 



50.5 51.0 51.5 52.0 52.5 

log(£'^et) erg 


Fig. 6.— Plot of the volumetric rate limit on choked GRBs in the nearby universe. The bin 
for canonical values of the RMW/AB emission model is marked by the star. The dashed 


line contour gives the rate of core-collapse supernovae within 10 Mpc as measured by Kistler 


et ah (2011). The dot-dashed line is the volumetric rate extracted from a large survey of 


SNe in the local universe (Leaman et ah 2011). 


The volumetric rate limit for a range of values of Ej and T;, is plotted in Figure]^ 
Two separate measurements of the nearby CC SNe are plotted as well to provide context 
to the calculated rate limits. Choked GRB events harboring particularly energetic jet 
parameters should be visible to the search method. However, if one compares the limits for 
the canonical RMW/AB model parameter values (F;, = 3, Ej = 10^^'^ erg) to the CC SNe 


















rates, it is clear that the current search method is not very sensitive to large regions of the 
model parameter space. However, the sensitivity of this search can be improved through 
rehnement of the event selection and analysis methods. Potential changes include greater 
signal retention through more efficient use of multi-variate machine learning cuts in the 
event selection process, the use of reconstruction methods optimized for sub-TeV muon 
tracks, and more accurate modeling of event angular error distribution. 


5. Conclusions 

The described search examined a newly developed data set consisting of 30-300 GeV 
muon neutrinos. No evidence for transient astrophysical neutrino sources was found in the 
data, leading to the construction of upper limits on the neutrino fluence of potential sources 
within the observation period. In particular, we examine the derived limit in the context 
of neutrino emission from choked GRBs. Although this search in its current conhguration 
is only sensitive to particularly energetic or nearby choked GRBs, the sensitivity of this 
method will improve as the event selection and search techniques are further optimized for 
muon neutrino events at sub-TeV energies. Gontinued development of this event selection 
will complement the current mature IceGube analyses at higher energies, leading to an 
overall enhancement of the detector’s sensitivity to transient sources. 
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